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Analysis  of  Numerical  Results 

We  begin  this  pan  by  desoilwg  a  uBijor  (SfTeience  between  circular  resonance 

_ 

orinfs  (geo^doonoas,  OPS),  and  an  elliptic  resonance  ort^  (MOLNIYA). 

As  described  in  ouf  previous  rqms,  the  disturinng  ftmcdcm  in  terms  of  orbital 
dements  s  giveii 

**S  X  X  x^« 

r=2  iHso  f =0  o=-o» 


where 


[(Z  -  2p)ffl + (Z  -  2p + q)M + m(Q  ~ 

HexeFhyrO')  and  0|p^(e}  are  the  inclinarion  and  eccentricity  functions  and  are 
explicitly  by  Raula  (l96d}«  Furthennoic.  $ = 0^  where  is  the  rotation  rate  of  the 

earth.  Since  10  the  first  approximation  M  =  ni,  fit »  constant  and  12  «  constant ,  die 
frequendes  are  given  by 


and  therefore  resonance  (/  »  0)  is  defined  by 


77t 


l-lp  +  q 


This  leads  to  the 


conclusion  diat  for  one-to-one  resonance  —  » I  any  combination  of  I  and  m  represems 


resonance.  For  two-ro-one  resonance 


(t-’) 


only  even  order  harmonics 


(«  =  2A6,...)  of  artttniy  degree  represent restmance. 


— 2*» 
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Older 

^*1 

degree 

order 

^*2 

degree 

m*2 

/=2,3.,.. 

m*2 

1*2,3,... 

Ma>4 

7*4,5.... 

m«3 

w*6 

7*S.7.... 

Hie  above  haimooies  lepies^t  zesonance  for  any  eccentria  edit  such  as  MOLNIYA  ozbit 
wth  —  =  2.  However,  for  circular  oxbits  (f  ■  0)  ooly  a  subset  of  the  above  hanuonics 

tepresrats  resonance.  These  subsets  me  dcGned  by  —  =  In  particular  the  subset 

no,  I’-lp 


t(x  iseosynchrsnoas  orHis 


"(i“} 


and  GPS  orbits  f- 


given  by 


Older 

degree 

order 

degree 

m*! 

7“lf3,5,7,.  ♦. 

?ft=2 

/=3,5,7.9.. 

«=2 

7*2,4,d,8,.*. 

m=4 

/*4,6,8,... 

»w=3 

7=3,5,7,9,... 

J7I=6 

7=7,9,,.. 

The  numerical  results  in  Ais  repoK  (see  Appendix)  fit  well  Ac  above  description.  Tbe 
average  effect  of  each  individual  harmonic  on  Ae  position  can  be  estimoied  from  the 
difference  of  two  coosecudve  PMS  values.  The  results  for  geosynchronous  and 
GPSotbite  are  pieseuted  in  Table  I  and  Table  II. 
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33lP,avgta^  fffrrn  nf  ,nrfi,>rduai 

20.^  DOSitinn  f/s.- .  ^  . 

1^2 

3- 

m 

1 

5 

6 

1  0.000 

0.792 

0.000 

0.002 

0.000 

2  7.973 

0.000 

0.092 

0.000 

0.000 

3 

0.278 

0.002 

•» 

0.013 

0.000 

4  _ _ 

— 

0.181 

0.000 

0.001 

S  * 

- - 

0.O55 

0.000 

6 - 

- - 

0.001 

Table  n 


Bmtsiiarfmi  ftf  MYirtmT  o„  h, — f- ,  ir  - 1  iTi|i 


1  ->  2 

T 

3 

4 

5 

6 

7 

1  0,000 

0.026 

0.000 

0.000 

0.000 

0.000 

2  0.822 

0.043 

0.007 

0.012 

0.000 

0.000 

3 - 

0.045 

0.015 

0.001 

o,ooo’ 

0.000 

4  ■ 

0.18S 

0.005 

0.003 

0.000 

5 - 

— 

— 

0.015 

0.003 

0.000 

6  —  , 

- 

0,000 

0.007 

7 - 

— 

— 

— 

0.000 
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Representation  of  Geopotential  with  Point  Masses 


Prepared  for  the  USAF  Office  of  Sdentific  Research 

by 

V.  Szdidicly,  K.  Zare,  and  S.  Chesley 


Concept 

The  standard  formulation  of  the  geopotendal  is  the  spherical  harmonic  expansion 
usixtg  Legendre  polynomials.  A  proposed  alternate  formularion  consists  of  representmg  the 
earth  by  a  large  number  of  poini  masses  positioned  on  a  fixed  grid  within  the  geotd.  We 
have  tiiree  primary  reasons  to  pursue  this  alternate  method: 

1 .  The  calculation  of  gravitational  acceleration  with  the  alternate  formulation  is  a  ample 
computation  using  Newton’s  law  of  gravity  for  eadi  mass.  The  same  calculation  with 
the  standard  Simulation  requires  the  recursive  computation  of  the  Legendre 
polynomials.  This  strongly  suggests  that  the  alternate  method  may  rignificantly  reduce 
the  computation  time  needed  to  provide  a  given  accura^.  (Or,  eqinvaleoily,  provide 
greater  accuracy  with  the  same  computational  eiqiense.)  Additionally,  the  point  mass 
method  is  ideally  suited  for  parallel  computation,  while  the  standard  formulation 
cannot  be  computed  in  parallel. 

2.  The  standard  formulation  has  a  slow  rate  of  convergence  due  to  resonant  terms  in  the 
harmonic  expansion.  This  leads  to  a  serious  problon  for  highly  eccentric  resonant  .  .  ■ 
orbits,  such  as  Molnya  orbits.  We  expect  H^er  rates  of  convergence  as  well  as  a  ' 
generally  monotonic  conv^gence  for  a  mass  point  model.  Hence  this  alternate  method 
should  prove  more  suitable  for  the  challenging  case  of  resonant  orbits. 

The  alternate  method  is  also  appealing  because  it  is  much  more  phy^cally  innutive.  An 
accurate  point  mass  model  would  by  its  very  nature  provide  the  three-dimensional 
density  distribution  of  the  earth.  Such  geophy^cal  information  may  be  used  to  improve 
the  modern  earth  models  and  it  would  be  of  interest  to  spedalists  outside  the  field  of 
orbital  mechanics. 
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1996  Accomplishments  previous  WortO 

l^elnitial  otgective  was  a  proof  of  the  concept,  v^di  was  accomplished  with  tiie 
follo\rag  steps.  ■ 

1 .  Assume  a  potential  model  as  the  *Hnith”  potential,  e.g.,  a  simple  potential  with 
1000  mass  points. 

2.  Generate  “tnith"  satellite  oihhs  from  the  potential  of  5t^  1 . 

3 .  Uaing  the  orbits  from  Step  2,  extract  **ob8ervations”.  These  observations  of  the 
^’truth’*  model  may  be  used  to  calculate  a  best  alternate  potential  model  in  a 
least  square  sense.  The  alternate  model  is  a  point  mass  model  with  100, 200, 
500,  etc.  total  mass  points. 

4.  Compare  the  orbits  generated  from  the  modd  potential  with  the  orbits 
generated  from  the  “mitb**  potential. 

As  was  reported  in  September  1 996,  tiie  proof  of  concept  phase  did  validate  or  refine 
most  of  the  basic  ideas  and  suggest^  a  clearer  future  course  of  action,  jij 

It  should  be  emphatozed  that  the  bulk  of  activity  in  the  first  yearwas  spent 
dtscussing  possible  stretches  and  implementing  these  strategies  on  a  computer.  The 
writing  and  testiiig  of  software  have  required  a  substantial  up-front  investment  of  lime, 
and  we  now  have  computa-  implem^itations  for  the  following  tasks. 

•  Create  a  point  mass  modd  with  a  prescribed  number  of  points.  The  modd 
consists  of  a  series  of  concentric  shells  of  pomts,  each  point  being  associated 
with  an  approximately  equal  volume.  The  mass  of  eadi  point  is  calculated 
based  on  a  con-unifotm  density  distribution  for  the  earth.  ' 

•  Generate  simulated  **truth"  observations  from  a  point  mass  or  spherical 
harmonic  potential.  These  can  be  obtained  for  a  number  of  different  satellites  in 
various  orbit  classes.  The  observations  may  come  in  varying  time  intervals  and 
can  have  random  noise  added  to  the  position. 

•  Estimate  the  potential  model  based  on  selected  observations,  hr  practice  one 
must  make  an  initial  guess  and  use  very  short  orbit  arcs  to  refine  that  guess. 
Due  to  the  poor  nature  of  die  initial  guess  a  long  orbit  arc  goes  beyond  the 
linear  region  and  the  method  fails.  As  the  guesses  become  more  and  more 
refined  the  orbit  arcs  used  as  observations  may  become  longer  and  longer,  tiius 
inqnoving  the  accuracy  of  the  estimate. 
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f997  Objectives  (Future  Worf^  . ,  . 

The  uhinate  ggal  of  this  research  is  to  obtain  a  working  mass  point  geopocential 
model  which  can  be  used  to  accuratdy  compute  oibhs  ofearth-bibiting  satellites.  There. 
are  two  preliminaiy  issues  that  need  to  be  considered  before  this  task  can  be  approached: 

•  The  computational  cost  of  estimadng  the  appropriate  mass  values  d^ends 
strongly  on  the  number  of  point  masses  used  in  the  model.  Given  a  deasity 
distribution  over  die  volume  of  die  earth,  how  many  mass  pomts  are  needed  to 
provide  a  given  levd  of  position  accurai!^  in  an  oibhal  simulation?  On  a  more 
prscdcal  level,  how  many  mass  points  are  needed  to  match  the  level  of 
aecuraqr  provided  by  a  given  tnmeation  of  the  standard  fbimulation,  e.g.,  a 
12x12  eq>ansioa?,71us  will  depend  on  the  complexly  of  the 
distribution,  and  to  a  lesser  extent  on  the  dass  of  orbit  being  analyzed. 

•  The  effiden^  of  the  estimation  process  is^  of  course,  influenced  by  the 
observational  data  provided  to  routme.  The  choice  of  the  number  of 
different  satellites  observed,  the  number  of  observations  for  each  satellite,  and 
the  time  interval  between  observations  (length  of  oibit  arc)  all  affect  the 
expense  of  the  estimation.  In  fact  a  poor  choice  can  lead  to  no  convergence.  Is 
there  a  reliable  way  of  determining  ^e  opthnal  choice  of  these  parameters 
before  die  estimation  has  begun?  If  not,  is  there  a  way  of  avoiding  tiie 
deddedlynon-optimal  choices? 

Once  these  questions  have  been  addressed  we  will  be  ready  to  use  precise  real-world  oibh 
ephemerides  to  calculate  a  working  geopotential  model.  Such  a  calculation  will  likely  ‘  ■ 
require  many  mass  points  udiieh  will  dictate  tiie  use  of  a  supocomputer.  One  riiould  note 
that  tiie  time  consuming  estimation  process  does  not  need  to  be  done  as  a  part  of  routme 
orbit  determination.  Once  a  "best”  point  mass  model  has  been  determined,  it  would  be  . .  - 
used  instead  of  the  standard  formulation  to  streamline  any  future  orbit  prediction  work. 
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